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ABSTRACT 

We report the discovery of a 6 1- Jupiter- mass brown dwarf, which transits its F8V host star, WASP- 
30, every 4.16 days. From a range of age indicators we estimate the system age to be 1-2 Gyr. We 
derive a radius (0.89 ± 0.02 i?j U p) for the companion that is consistent with that predicted (0.914 
Rjup) by a model of a 1-Gyr-old, non-irradiated brown dwarf with a dusty atmosphere. The location 
of WASP-30b in the m i nimum of the mass-radius relation is consistent with the quantitative prediction 
of IChabrier fc Bara ffc (2000]), thus confirming the theory. 

Subject headings: binaries: eclipsing — stars: low-mass, brown dwarfs — stars: individual (WASP-30) 



1. INTRODUCTION 

A brown dwarf (BD) is traditionally defined as an 
object with a mass above the deuteriu m-burning limit 
(13 M Jup ; e.g. IChabrier etafl l2000a| ) and below the 
hydrog en-burning limit (0.07 M ; e.g. IChabrier et al.l 
l2000bf ). However, an alternative suggestion is that 
the manner in which an object forms should determine 
whether it is a planet or a BD. Thus, if an object formed 
by core accretion of dust and ices in a protoplanetary 
disc then it would be a planet, and if it formed by gravo- 
turbulent collapse of a molecular cloud, as do stars, then 
it would be a BD. 

Studies such as the iCaballero et al.l (j2007) observa- 
tions of a young open cluster core find a continuous 
mass function down to ~6 -Mj up , indicating that the 
star formation mechanism can produce objects with plan- 
etar y masses. This is suppor t ed by theoretical stud- 
ies (Padoan fc NordluncH [2001 iHennebelle fc Chabrierl 
12008ft which suggest that gravoturbulent fragmentation 
of molecular clouds produces stars and BDs down to 
a few Jupiter masses in numbers comparable to the 
observationally-determined distribution. In contrast, 
when taking into account planetary migration through 
the protoplanetary disc, the core accretion process might 
result in giant plan ets with masses of up to 10 M,t„ p 
(lAlibert et al.ll2005D. or even 25 M Jup (jMordasini et alj 
2008). Sahlmann et al.l (|2010l ) see evidence for a a bi- 
modal distribution in brown dwarf masses, with the less- 
massive group presumably representing the high-mass 



dra@astro.keele.ac.uk 

1 Astrophysics Group, Keele University, Staffordshire, ST5 
5BG, UK 

2 SUPA, School of Physics and Astronomy, University of St. 
Andrews, North Haugh, Fife, KY16 9SS, UK 

3 Observatoire de Geneve, Universite de Geneve, 51 Chemin 
des Maillettes, 1290 Sauverny, Switzerland 

4 Astrophysics Research Centre, School of Mathematics & 
Physics, Queen's University, University Road, Belfast, BT7 INN, 
UK 

5 Department of Physics and Astronomy, University of Leices- 
ter Leicester, LEI 7RH, UK 

° Institut d'Astrophysique et de Geophysique, Universite de 
Liege, Allee du 6 Aout, 17, Bat. B5C, Liege 1, Belgium 

7 Las Cumbres Observatory, 6740 Cortona Dr. Suite 102, 
Santa Barbara, CA 93117, USA 



tail of the planetary distribution. 

An accurate, precise measurement of an object's radius 
is therefore required to probe for the existence of a core 
and thus discriminate between the two formation mech- 
anisms. For example, the radius of the 8-Mj up body, 
HAT-P-2b, is consistent with an irradiated planet in- 
corporating a 340- Earth-mass core, but is smaller than 
if it were c oreless dLeconte et al.l 120091). The 22-M Jup 
CoRoT-3b (|Deleuil et al.ll2008f ) is sufficiently massive to 
qualify as a BD under the traditional definition, but the 
radius of this object is uncertain at the 7% level. This 
is higher than the 3% required to discriminate between 
the absence or presence of a core a nd thus det e rmine 
how it formed (jLeconte et al.1 120091 ) . Ilrwin et al.l (j2010f ) 
found a ~30-M Jup BD, NLTT 41135C, which transits 
one member of an M-dwarf binary system. However, as 
the transits are grazing, it is not currently possible to 
accurately measure its radius. 

There is less ambiguity around the upper end of the 
BD mass regime: if a body is sufficiently massive to fuse 
hydrogen then it is a star, otherwise it is a BD. High- 
mass BDs with precise radius measurements are useful 
for testing BD evolution models, as it is in the high-mass 
regime that m odels predict the gre a test changes in radius 
with age (e.g. [Baraffc ct al. 2003). Stass un et al.l ()2006[ ) 
discovered a BD eclipsing binary system in the Orion 
Nebula star-forming region, with masses of 57 ± 5 Mj up 
and 36 ± 3 M Jup . With very large radii of 0.699 ± 0.034 
Rq and 0.511±0.026 i?©, it seems that these objects are 
in the earliest stages of gravitational contr action. Similar 
to th e NLTT 41135 system, LHS 6343 C (pohnson et al.l 
I2010L Johnson 2010, private communication) is a 63- 
Mj up BD that transits one member of an M-dwarf bi- 
nary system. In this case, the transits are full and so 
the radius (0.825 ± 0.023 i? j UP ) of this object is precisely 
determined. CoRoT-15b (Bouchv et al.l 120101) is a 63- 



Mj up -mass BD in a 3-d orbit around an F7V star. Due 
to the faintness of the host star (V ~ 16), the BD radius 
(1.12j. ; 15 iijup) is not yet well determined. 

To test and refine models of BD formation and evo- 
lution, a population of well-characterised objects is re- 
quired. In this letter, we present the discovery of WASP- 
30b, a 61-Mj up brown dwarf that transits its moderately- 
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TABLE 1 

Radial velocity measurements of WASP-30 
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bright host star. 



2. OBSERVATIONS 

WASP-30 is a V = 11.9, F8V star located in Aquar- 
ius, on the border with C etus. A transit search 
(|Collier Cameron et al.ll200ll of WASP-South data from 
2008 July to November found a strong, 4.16-d periodic- 
ity. Furthe r observations in 20 09 with both WASP in- 
struments (|Pollacco et al.l 120061) led to a total of 17612 
usable photometric measurements (Figure [1} . 

Using the CORALIE spectro graph mounted on 
the 1.2-m Euler-S wiss telescope (jBaranne et al.l 119961 : 
lOueloz et ai1l2000D . we obtained 16 spectra of WASP-30 
in 2009. Radial velocity (RV) mea surements were com- 
puted by weighte d cross-correlation (jBaranne et al.lll996l : 
Pepc ct al. 2005) with a numerical G2-spectral template. 
RV variations were detected with the same period found 
from the WASP photometry and with a semi-amplitude 
of 6.6 km s _1 , consistent with a sub-stellar-mass com- 
panion. The RV measurements are listed in Tableland 
are plotted in Figure [TJ 

To test the hypothesis that the RV variations are due 
to spectral line distortions caused by a blended eclipsing 
binary or starspots, we performed a line-bisector analysis 
(|Queioz et al.l l200F) of the CORALIE cross-correlation 
functions. The lack of correlation between bisector span 
and RV (Figure [lj supports our conclusion that the peri- 
odic dimming of WASP-30 and its RV variations are due 
to the sub-stellar orbiting body, WASP-30b. 

To refine the system parameters, we obtained high 
signal-to-noise transit photometry through a Gunn r fil- 
ter with the Euler-Swiss telescope on 2010 Aug 01 (Fig- 
ure [2 Table [2]) . The data were affected by light cloud 
and a guiding issue caused by the close proximity of the 
bright Moon (69% illumination, 26° from target). 

3. STELLAR PARAMETERS 

The 16 CORALIE spectra of WASP-30 were co-added 
to produce a single spectrum with a typical S/N of 
around 70:1 . The analysis was p erformed using the meth- 
ods given in lGillon et al.l (|2009ft . The H a line was used to 
determine the effective temperature (T G ft), while the Na 
I D and Mg I b lines were used as surface gravity (log g*) 
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Fig. 1. — Top panel: WASP discovery light curve, folded on 
the ephemeris of Table \3\ Second panel: Binned (A</> = 0.0025) 
WASP data around the transit. Third panel: Transit light curve 
obtained with Euler. Fourth panel: relative RV measurements of 
WASP-30 measured by CORALIE. Fifth panel: residuals about 
the model RV solution. Bottom panel: bisector spans. 
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TABLE 2 




EULER PHOTOMETRY OF WASP-30 
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5409.693659 
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This table is published in its entirety in the electronic edition of 
the Astrophysical Journal. A portion is shown here for guidance 
regarding its form and content. 

diagnostics. The parameters obtained from the analy- 
sis are given in the top panel of Table [3J The elemen- 
tal abundances were determined from equivalent width 
measurements of several clean and unblended lines. A 
value for microturbul ence (&) was d etermined from Fe I 
using the method of iMagainl (|1984f ) . The quoted error 
estimates include that given by the uncertainties in T g, 
logg* and £t, as well as the scatter due to measurement 
and atomic data uncertainties. Our quoted lithium abun- 
dance takes accou nt of non local therm odynamic equilib- 
rium corrections ( [Carlsson et al.|[l994f) . with a value of 
Au = 2.95 resulting when neglecting them. 

The projected stellar rotation velocity (u sini) was de- 
termined by fitting the profiles of several unblended Fe I 
lines. A value for macroturbulence (f m ac) of 4.7 ± 0.3 
km s _1 was assumed, based on the tabulation by iGravl 
(|2008D . and an instrumental FWHM of 0.11 ± 0.01 A, 
determined from the telluric lines around 6300A. A best- 
fitting value of vs'mi = 14.2 ± 1.1 km s _1 was obtained. 

4. SYSTEM PARAMETERS 

The WASP and Euler photometry were combined 
with the CORALIE RV measurements in a simulta- 
neous Markov-chain Monte Carlo (MCMC) an alysis 
(|Collier Cameron et al.l [20071: iPollacco et al.| [2008). Our 
proposal parameters are: T c , P AF, T14, b, K\, T e s 
[Fe/H], y/ecosa; and \ fe sin w (|Collier Cameron et al.l 
[20071 lEnoch et al.|[2010h . Here T c is the epoch of mid- 
transit, P is the orbital period, AF = Rp/Rl is the frac- 
tional flux-deficit that would be observed during transit 
in the absence of limb-darkening, T14 is the total transit 
duration (from first to fourth contact), b is the impact 
parameter of the BD's path across the stellar disc, K\ 
is the semi- amplitude of the stellar reflex velocity, T e g 
is the stellar effective temperature, [Fe/H] is the stellar 
metallicity, e is the orbital eccentricity and w is the ar- 
gume nt of periast ron. 

As iFordl ([2001 notes, it is convenient to use ecosw 
and esincj as MCMC jump parameters, because these 
two quantities are nearly orthogonal and their joint prob- 
ability density function is well-behaved when the eccen- 
tricity is small and u) is highly uncertain. Ford cautions, 
however, that the use of ecosw and esinw as jump pa- 
rameters implicitly imposes a prior on the eccentricity 
that increases linearly with e. Instead we use -^/ecosw 
and y/e&inuj as jump parameters, which restores a uni- 
form prior on e. 

At each step in the MCMC procedure, each proposal 
parameter is perturbed from its previous value by a 





""""^^^ '■ 

0.2 Gyr 


's 

s' 


1 .0 Gyr 




2.3 Gyr 




— - 5.0 Gyr 

I 





0.2 0.4 

J-H 

1.6 



1.4 



Q. 

^ 1.2 

CO 



0.8 



0.6 

Mass / M Jup 

Fig. 2. — Top panel: Colour-magnitude diagram for WASP- 
30 and a nearby star with which it appears to be com oving. 
Isochrones for the ages shown are from Mari go et all (2008 ) . Bot- 
tom panel: Mass-radius diagram showing objects in the mass 
range 5-110 Mj up with precise radius measurements. The data 
are from http://www.inscience.ch/transits/ and, in the case of 
LHS 6343 (J, from Johnson (2010, private communication). The 
theoretical isochrones with ages 0.1, 0.5, 1, 5 and 10 Gyr are from 
the DUSTY00 and COND03 models. Those that reach down to the 
minimum plotted mass are the COND03 isochrones. The vertical 
dashed lines depict the approximate theoretical minimum masses 
for deuterium burning (13 Mj up ), lithium burning (0.055 Mq), 
and hydrogren burning (0.07 Mq; e.g. IChabrier et al.|[2u00al lb|l. 



small, random amount. From the proposal parameters, 
model light and RV curves are generated and x 2 is cal- 
culated from their comparison with the data. A step 
is accepted if x 2 ( our merit function) is lower than for 
the previous step, and a step with higher \ 2 is accepted 
with probability exp(— A% 2 /2). In this way, the param- 
eter space around the optimum solution is thoroughly 
explored. The value and uncertainty for each parameter 
are respectively taken as the median and central 68.3% 
confidence interval of the parameter's marginalised pos- 
terior probability distribution. 

From the proposal parameters, we calculate the mass 
M, radius R, density p, and surface gravity log g of the 
star (which we denote with subscript *) and the planet 
(which we denote with subscript P). At each step, the 
stellar density is measured from the transit light curve 
(jSeaeer fe Mallen-Ornelai |2003[ ). This is inpu t in t o 
the empirical mass calibration of iTorres et al.l (2010), 
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as modified by lEnoch et al.l (|2010T ) . to obtain an esti- 
mate of the stellar mass. We also calculate the equilib- 
rium temperature of the planet T eq i, assuming it to be 
a black-body with efficient redistribution of energy from 
the planet's day-side to its night-side, the transit ingress 
and egress durations, T\i and T34, and the orbital semi- 
major axis a. 

With eccentricity flo ating, we find e = 0.002 1^°;°^. 
Applying the 'F-test' of lLucv fc Sweenevl (|1971h . we find 
a 70% probability that the fitted eccentricity could have 
arisen by chance if the the underlying orbit is in fact 
circular. As such, we impose a circular orbit, but we note 
that doing so has no significant effect on the solution in 
this case. 

Without exquisite photometry, our implcntation of 
MCMC tends to bias the impact parameter, and thus 
i?» and i?p, to higher values. This is because, with low 
signal-to-noise photometry, the transit ingress and egress 
durations are uncertain, and symmetric uncertainties in 
those translate into asymmetric uncertainties in b and 
thus on We therefore place a main sequence (MS) 
prior on the star, which is reasonable given the star's ap- 
parent age (Section [5]). With the MS prior, a Bayesian 
penalty ensures that, in accepted MCMC steps, the val- 
ues of stellar radius are consistent with the values of stel- 
lar m ass for a main-sequence star (|Col licr Ca meron et al.l 
l2007h . 

The median parameter values and their 1-a uncertain- 
ties from our MCMC analysis are presented in the mid- 
dle panel of Table [3] The corresponding transit light 
curves and RV curve are shown in Figure [TJ When not 
imposing a MS prior, the best-fitting values we obtain 
are: b = 0.24^, = 1.337±° ^ R Q , and R P = 
0.925_q -Rjup- 

5. SYSTEM AGE AND COMPANION RADIUS 

The high lithium abundance (A u = 2.87±0.I0) found 
in WASP-30 implies an age most likely between that of 
open clusters such as a Per (50 Myr; Au — 2.97 ± 0.13) 
and the Hyades (600 Myr; A Li = 2 .77±0.21), and almost 
certainly younger than 2-Gyr- old open clusters such a s 
NGC 752 (Au = 2.65 ±0.13; ISestito fc Randichl 120051) . 

Assuming aligned stellar-spin and planetary-orbit axes, 
the measured vsini of WASP-30 and its derived stel- 
lar radius indicate a rotational period of P ro t = 4.6 ± 
0.4 days. After removing the transits from the WASP 
light curves, we searched them for evidence of rotational 
modulation. Though there are periodogram peaks at pe- 
riods of 4.1 d, 4.3 d and 4.7 d, the signal amplitudes 
are small. We thus conclude that there is no evidence of 
rotational modulation in the WASP-30 light curves, com- 
mensurate with expectations based on the star's spectral 
type. The 4.6 d stellar rotation period is very close to 
the companion's orbital period, suggesting that the two 
may be synchronized, t hus preventin g a gyrochronologi- 
cal age determ inatio n (Barnes 2007). The synchroniza- 
tion timescale (|Zahnlll977t ) for the star is 0.23±0.02 Gyr. 

We searched within 15' of the sky position of WASP- 
30 for common proper motion stars. The V = 13.6 
star USNO-B1.0 0800-0674908 is 13f09 away and ap- 
pears to be comoving with WASP-30 (Table |31 bot- 
tom panel). Using 2MASS photometry, we constructed 
a colour-magnitude diagram. A distance modulus of 



TABLE 3 
System parameters for WASP-30 



Stellar parameters from spectroscopic analysis 
R.A. = 23 h 53 m 38.03 s , Dec. = -10° 07' 05. l" (J2000) 
TYC 5834-95-1, 2MASS 23533805-1007049 

6100 ± 100 
4.3 ± 0.1 



T eS (K) 
log g, 
it (km s" 1 ) 
v sin i (km s^ 1 ) 
[Fe/H] 
[Si/H] 
[Ca/H] 
[Ti/H] 
[Ni/H] 



1.1 ± 0.1 
14.2 ± 1.1 
-0.08 ± 0.10 
+0.04 ± 0.13 
+0.10 ± 0.14 
-0.01 ± 0.14 
-0.08 ± 0.13 
2.87 ± 0.10 



Parameters from MCMC analysis 



P(d) 

T c (HJD) 

T14 (d) 

T12 = T 34 (d) 

AF = Rl/Rl 

b 

in 

K ± (km s- 1 ) 

a (AU) 

e 

7 (km s _1 ) 

Af* (M ) 
R, (P ) 
log9* (cgs) 

p* (pe) 

T efI (K) 

[Fe/H] 

A/ P (Mjup) 
Rp (Rju P ) 
loggp (cgs) 
Pp (pj) 
T cq i (K) 



4.156736 + 0.000013 
2455334.98479 ± 0.00076 
0.1595 + 0.0017 
0.01060 ± 0.00025 
0.00498 ± 0.00017 



0.066 



+0.072 
-0.044 



89.57 



+0.28 
■0.47 



6.627 + 0.015 
0.05325 ± 0.00039 
(adopted) 
7.9177 + 0.0099 

1.166 + 0.026 
1.295 + 0.019 
4.280 + 0.010 
0.537 + 0.019 

6201 ± 97 
-0.03 + 0.10 

60.96 ± 0.89 
0.889 ± 0.021 
5.247 + 0.019 
86.8 + 5.7 
1474 ± 25 



Proper motions of WASP-30 and its neighbour 



Star 

WASP-30 
neighbour 



UCAC3 PPXML 
Ha (mas/yr) (mas/yr) /j q (mas/yr) fj, s (mas/yr) 



-20.4 + 2.2 
-23.0 + 7.2 



-9.0 + 2.4 
-8.1 + 4.4 



-22.9 + 2.2 
-18.9 + 3.9 



-9.0 + 2.4 
-9.5 ± 3.9 



8.50±0.05 (500±10 pc) was required to place the comov- 
ing star on the main-sequence and suggests that WASP- 
30 is ~2.3-Gyr old (Figure [U upper panel). However, 
the apparent magnitude and spectral type of WASP-30 
suggest a smaller distance modulus of 7.9 ±0.2 (366 ± 77 
pc). As the comoving star may be a mere line-of-sight 
neighbour, this age determination should be treated with 
caution. 

In the lower panel of Figure [U WASP-30b is plot- 
ted in a mass-radius diagram together with isochrones 
of isolated BDs from mode l s with dusty atmospheres 
(DUSTY00, iChabrier et all l2000b[> and models wit h 
dust- free atmospheres (COND03. iBaraffe et al.l 12003). 
By an age of 1 Gyr, an isolated BD with the mass of 
WASP-30b is expected to have cooled to T cql ~ 1700 
K, and the transition from dusty L-dwarfs to dust-free 
T- dwarfs is expected t o take place at T eq \ = 1300-1700 
K (|Baraffe et al.l 12003). Considering this and the fact 



WASP-30b: a 61 M Jup brown dwarf 



5 



that WASP-30b is highly irr adiated, it is likely that the 
dusty atmosphere models of iChabrier et al.l (|2000bD are 
more representative. Depending on the transition tem- 
perature, WASP-30b may remain a dusty L-dwarf for the 
lifetime of its host star, or it may at some point transition 
to a dust-free T-dwarf. 

WASP-30b has a mass of 0.05819 ± 0.00084 M© and 
a radius of 0.0914 ± 0.0022 R G . The DUSTY00 mod- 
els predict the radius of an isolated, 0.06-Af Q BD to be 
0.149, 0.104, 0.094 and 0.084 i? at ages of, respectively, 
0.1, 0.5, 1 and 5 Gyr. From a simple linear interpola- 
tion of the DUSTY00 model values, the measured ra- 
dius of WASP-30b and its 1-rr uncertainty suggests its 
age is 2 ± 1 Gyr. The measured radius is inconsistent 
with the DUSTY00 model for a 0.5-Gyr BD at the 5.7- 
cr level and inconsistent with a 5-Gyr BD at the 3.4-ct 
level, but consistent for a 1-Gyr BD at the 1.2-er level. 
This agreement would be slightly be tter if the DU S TY00 
models took account of irradiation. iBaraffe et al.l (|2003[ ) 
find that the irradiation of a dust-free atmosphere, at the 
level of irradiation experienced by WASP-30b, results in 
radii larger by 10% for a 1-Mj up planet, and larger by 
7% for a 10-Mj up planet. However, they note that the 
evolution of dust-free atmospheres are more affected by 
irradiation than are dusty atmospheres, and WASP-30b 
is considerably more massive. 

WASP-30's lithium abundance favours a young system 
age of 50-600 Myr, though lithium is a poor age indicator 
for an F8 star, and an age of up to 2 Gyr is not ruled out. 
The apparent rotational synchronisation of the host star 
places a lower limit of 200 Myr (though the system may 
have been born synchronised) and a possible companion 
star suggests an older age of 2.3 Gyr. Given the BD's 
measured radius, the DUSTY00 BD model indicates an 
age of 2± 1 Gyr. Taken together, we suggest that an age 
of 1-2 Gyr is most likely. 



the MRR at high BD masses (see Section 3.1 of that 
paper). The location of WASP-30b in the MRR min- 
imum is consistent with the quantitative prediction of 
IChabrier" fc Baraffc (2000), thus confirming the theory. 

Thus far, we know of two other high-mass brown 
dwarfs that transi t stars: CoRoT-15b (63 -Mj up ; 
Bouchvet all f2010h and LHS 6343 C (63 Mj up ; 
Johnson et al.l 120101 Johnson 2010, private communica- 
tion ). The radius of CoRoT-15b (1.121^° i? Jup ) is 
uncertain and the age of the system is currently un- 
constrained. The faintness of the host star (V ~ 16) 
makes improving this situation difficult. LHS 6343 G 
was found to transit one member of an M-dwarf binary 
system using Kepl er photometry (KIC 10002261; e.g. 
IBorucki et aL|[2010[ ). It initially seemed that LHS 6343 C 
was larger than pred icted for a BD of its mass and age 
(j Johnson et al.l[20~Toh . However, after a re-evaluation of 
the treatment of the third light in the system, it seems 
t o be consistent ( J ohnson 2010, private communication). 

IChabrier et alj (|2000bD predict that BDs with M < 
0.05 M Q do not burn lithium, those with M > 0.06 M 
burn essentially all their lithium by an age of 0.5 Gyr, 
and those with an intermediate mass (M = 0.055 Mq) 
burn half their lithium by 0.5 Gyr, two-thirds by 1 Gyr, 
and three-quarters by 5 Gyr. With a mass of 0.0582 ± 
0.0008 and an age of 1-2 Gyr, WASP-30b is likely to have 
burned most of its supply of lithium. 

WASP-30b has the second smallest companion-to-star 
size ratio (AF = R*/Rl = 0.0050) of all sub-stellar 
bodies so far discovered by ground-based transit surveys. 
The star in the system with the smaller size ratio, HAT- 
P-ll (V = 9.6; AF = 0.0033: iBakos et al.ll201ff) . is 8 
times brighter than WASP-30. As it is far easier to find 
such an object around a smaller, cooler star, the discov- 
ery of WASP-30b suggests that high-mass, sub-stellar 
objects in short orbits around cooler stars are rare. 



6. DISCUSSION 

The discovery of WASP-30b heralds the first unam- 
biguous observational determination of the mass-radius 
relation (MRR) in the BD regime, and so we have added 
BDs to white dwarfs and neutron stars in the list of 
q uantum-dominated objects with radius determinations. 

IChabrier fc Baraffe! 1)20001 ) perfomed the first quanti- 
tative theoretical calculation of the MRR in the substel- 
lar and low-mass star domain, predicting a minimum in 
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the Issac Newton Group on La Palma. We are grate- 
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